Local density approximation for long-range or for short-range energy functionals? 
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Density functional methods were developed, in which the Coulomb electron-electron interaction is 
split into a long- and a short-range part. In such methods, one term is calculated using traditional 
density functional approximations, like the local density approximation. The present paper tries 
to shed some light upon the best way to do it by comparing the accuracy of the local density 
approximation with accurate results for the He atom. 



I. INTRODUCTION 

In recent years, there has been a growing interest in ap- 
proaches of density functional theory (DFT) yj based on 
a long-range/short-range decomposition of the Coulomb 
electron-electron interaction. The idea is to use different, 
appropriate approximations for the long-range and the 
short-range contributions to the exchange and/or corre- 
lation energy density functionals of the Kohn-Sham (KS) 
scheme Q. 

Various schemes combining a (semi)local short-range 
exchange energy functional approximation with an "ex- 
act" long-range exchange energy expression have been 
proposed (see, e.g., Refs. H 0, H II IS 0), which al- 
low to impose the exact Coulombic asymptotic behav- 
ior 1/r of the exchange interaction. This turned out 
to be important for charge transfer, van der Waals in- 
teractions, etc. But opposite approaches combining a 
(semi)local long-range exchange functional approxima- 
tion with an "exact" (or "hybrid") short-range exch ang e 
have also been used (see, e.g., Refs. H [H], 0, IT^) . 
which allow to introduce exact exchange in solid-state 
calculations without the computationally-demanding ex- 
act long-range contribution. There is also a physical 
reason for doing it: long-range correlations are not well 
treated by (semi)local density approximations. Treating 
only the exchange exactly destroys the balance of errors, 
which is important for metals, gaps, etc. Contemplating 
this two opposite approaches, one can ask: is it prefer- 
able to use (semi)local density functional approximations 
for the long-range or for the short-range contribution to 
the exchange energy? 

Nozieres and Pines EH have first used the idea of de- 
composing the correlation energy of the uniform elec- 
tron gas into long-range and short-range contributions 
to facilitate its calculation. In the context of DFT, a 
few schemes combining a (semi)local short-range density 
functional approximation with a long-range correlation 
energy calculated by other means have been proposed 
(see, e.g., Refs. El El E3 El El to deals with for ex- 



ample near-degeneracy or long-range van der Waals inter- 
actions. Indeed, it is well known that (semi)local density 
functional approximations are appropriate for the short- 
range contribution (see, e.g., Refs. l20t l2ll) . However, for 
a given decomposition of the Coulomb interaction, there 
are at least two possible definitions for a short-range cor- 
relation functional. It can be defined either as the differ- 
ence between the Coulombic correlation functional and a 
long-range correlation functional associated to the long- 
range part of interaction (as in Ref. Il7j) . or directly from 
the short-range part of the interaction (as in Ref. 1 1 6T> . 
Therefore, one can wonder: what is the preferable defi- 
nition for a short-range correlation functional as regards 
as the accuracy of (semi)local density functional approx- 
imations? 

This work sheds some light on these two questions. 
Taking the example of the He atom, we test the accu- 
racy of the local density approximation (LDA) to the 
long-range and short-range contributions to the exchange 
and correlation energies for a given decomposition of the 
interaction. Atomic units (a.u.) are used throughout this 
paper. 



II. LONG-RANGE AND SHORT-RANGE 
DENSITY FUNCTIONALS 

We define in this section the long-range and short- 
range density functionals that we consider in this work. 

Let's define first a general exchange functional for an 
arbitrary electron-electron interaction w ee (r) 



E x [n;w ee ] = ($[n]|Wee|$M) - E H [n;w ee ], 



(1) 



where W ee = (1/2) JJ n 2 {ri, r 2 )w ee (ri 2 )dri<ir2 is the 
interaction operator, expressed with the pair-density 
operator n.2(ri,r2), 3>[n] is the KS determinant and 
E-n[n;w ec ] = (1/2) jj n(r 1 )n(r 2 )w ee (r 1 2)dr 1 dr2 is the 
Hartree functional. Likewise, the general correlation 
functional associated to w ee (r) writes 

E c [n;w ee ] = min(*|f +W ee \V) - (*[n]|f + W ee \$[n}), 



(2) 
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where the constrained-search formalism has been 
used. 
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We use the following decomposition of the Coulomb 



electron-electron interaction 



V) = l/r 



(3) 



where wj^(r) = cvi{(ir)/r is a long-range interaction 
and = erfc(/zr)/r is its short-range complement. 

The parameter (i which controls the range of the de- 
composition varies between and oo. For fi = 0, the 



long-range interaction vanishes, w\ 



lr,/j=0 



0, and the 



short-range interaction reduces to the Coulomb interac- 



tion wll'^ 



w. 



coul 



In the limit (i — > oo, the long- 



range interaction reduces to the Coulomb interaction 
yjir./i-xxj _ w coui^ an( j ^. ne s hort-range interaction van- 
ishes, < r e '^°° = 0. 

The Coulombic exchange functional £"™ ul [7i] = 
E x [n; w™ 1 ] is consequently decomposed as 



e: 



coul r 



E*>"[ n ] + E»«[ n ], 



(4) 



where E^^[n] = E x [n; w *^] is the long-range ex- 
change functional associated to the interaction toj,'/ , and 
E^'^[n] = E x [n; wf^] is the short-range exchange func- 
tional associated to the interaction wf^. 

One can also define a long-range correlation functional 
associated to the interaction w^, E^'^[n] = E c [n; toj, 1 '/], 
and a short-range correlation functional associated to the 
interaction wf^, E s J^[n] — E c [n; wf^] . However, be- 
cause E c [n;w ee ] is not linear with respect to w ee , we 
do not have for the Coulombic correlation functional 
^coul^j _ E c [n; wl° ul ] the same decomposition as in 
Eq. (jtj) 



E?*[n]?E*>»[n]+E?>»[n]. 



(5) 



Instead, one can write i?™ ul [n] exactly as 

Er l M = E^[n] + E s c ^[n] + E^[n], (6) 

defining the mixed long-range/short-range correlation 
functional i?' r_sr, ' 1 [n] which encompasses all the terms 
stemming from the non-linearity of E c [n;w ee ] with re- 
spect to w ee . By associating the mixed term E l J~ S1 '^ [n] 
to either the long- or the short-range part of the corre- 
lation energy, one obtains two possible long-range/short- 
range decompositions. The first one is 



e: 



coul r 



E^[n]+E^[n}, 



(7) 



defining the new short-range correlation functional 
Ef'^[n] = Ef^ l [n] + E^~ sr '^[n}. The second one is 



El 



(8) 



defining the new long-range correlation functional 
E^[n] = E^[n] + E l *- SI >»[n]. 

All the long-range functionals vanish for (i = 0, 

££, M =0 = £lr, M =0 = £lr, M =0 = Q> and reduce to tne 

Coulombic functionals for (i —> oo, E 1 *'^ 00 = E x oul and 
E h,^oc = E h^oo = ^coui Symmetrically, all the 



short-range functionals reduce to the Coulombic func- 
tional for (1 = 0, Ef'^ = E™ ul and Ef<^° = 
Ef>» =0 = E™ nl , and vanish for (i -> oo, E£>^°° = 
Ef'^°° = Ef'^°° = 0. The mixed long-range/short- 
range correlation functional vanishes for both (i — and 
(j, oo, £*-«-./*=° = El 1 - 51 '^ 00 = 0. 

At first sight, the physical meaning of the decomposi- 
tion of Eq. © may not be obvious. To get more insight 
on the nature of the terms E 1 *' 1 *, Ef-^ and ££ r ~ sr '' i , we 
show now the expressions of these functionals in second- 
order Gorling-Levy perturbation theory J23|. We start 
from the second-order Coulombic correlation energy 



£Coul,(2) = ^ 



V'COUl 
Hx 



E s 



(9) 



where Vg° ul 



J n(r)5E^" 1 1 'Sn(r)dr is the Coulombic 
Hartree-exchange potential operator expressed with den- 
sity operator n(r) and with the Coulombic Hartree- 
exchange functional E^ 1 [n] , $ and E s are the KS wave 
function and energy, <E>i and E s a are the excited KS eigen- 
functions and eigenvalues. Applying in Eq. the long- 
range/short-range decomposition on IF™ ul and V£° ul and 
expanding leads to the second-order expressions for Eg' 1 *, 
Ef* and E^- S1 ^: 

• The second-order long-range correlation energy 
writes 



i 



where V, 



E.„ — E s 



(10) 



J h(r)SE^ /Sn(r)dr is the long- 
range Hartree-exchange potential operator ex- 
pressed with the long-range Hartree-exchange func- 
tional E 1 ^ [n] defined in the same way as the long- 
range exchange functional. 

• Likewise, the second-order short-range correlation 
energy writes 



E s — E s a 



(11) 



where — J n(r)5E^ /5n(r)dr is the short- 

range Hartree-exchange potential operator ex- 
pressed with the short-range Hartree-exchange 
functional ^HaTW defined in the same way as the 
short-range exchange functional. 

• Finally, the second-order mixed long-range/short- 
range correlation functional is 

R lr-sr,/i,(2) . 

2 y (gjjfr^ ~ - V£f\*) _ (12) 

E„ — E„a 
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III. LOCAL DENSITY APPROXIMATION 

For a given electron-electron interaction, the local den- 
sity approximation to the previously introduced exchange 
and correlation functionals consists in locally transferring 
the corresponding energy of a uniform electron gas with 
the same interaction and with density equal to the local 
value of the inhomogeneous density. 

Let us first consider the short-range functionals. 
The short-range exchange LDA functional associated to 
Ef'^ln] writes 

CiTdaM = / ™(r)4XifWr))dr, (13) 

where s s ^'^ nil (n) ^ s ^ ne exchange energy per particle of 
a uniform electron gas with interaction vf^ 11 p3. |25| . 
Similarly, the short-range correlation LDA functional for 

Ef'^n] is 

^l"daN = / n(r)t^ ni Mr))dr, (14) 

where g^uf( n ) = £™uif ( n ) ~ 4?unif ( n ) is obtained as the 
difference of the Coulomb correlation energy, e^^(n), 
and of the long-range correlation energy of a uniform 
electron gas with interaction w 1 ^ The LDA for 

Ef'^n] writes 

K&M = J (15) 

where s s J'^ nii is the correlation energy per particle of a 
uniform electron gas with interaction w|g' M |2t| . 

The LDA corresponding to the long-range function- 
als are obtained by difference to the Coulombic case: 

^cTdaM ~ 4TdaM and = £ c co l u daM " 

^ldaW' 

IV. RESULTS FOR THE HE ATOM 

We consider the simple example of the He atom. We 
compare the accuracy of the LDA for the long-range or 
short-range energies as follows. For each value of mu 
(selected in [0, oo|_), we obtain an accurate value of the 
energy (see Refs. ll7l.l27l for details), we calculate the LDA 
error on the energy, and report in the plots the LDA error 
(ordinate) for each accurate energy (abscissa) . We choose 
this way of plotting in order to emphasize the importance 
of the error for a given energy value to be recovered by 
the approximation: a good functional would yield small 
errors even for large contributions of the energy. Further- 
more, this representation gives a unique scale for both 
long- and short-range expressions. The left end of the 
plots corresponds to vanishing functionals (i.e., /i = for 




0.2 0.4 0.6 0.8 1 

-EjJ (a.u.) 



FIG. 1: LDA errors on E* 1 '" (solid curve) and E%'" (dashed 
curve) with respect to — i?| r ' M and —E 1 *'* 1 respectively, for the 
He atom. 



long-range functionals and fj, — > oo for short-range func- 
tionals). The right end of the plots corresponds to the 
usual (Coulombic) LDA case. 

Let us consider first the LDA errors on E^^ and E^'^ 1 , 
represented in Fig. ^ For all values of the exchange en- 
ergy (except of course at the end points) , the LDA error 
on E^'^ is systematically smaller than the error on E 1 *'^ 
(the difference is of the order of 0.025 a.u. at an inter- 
mediate energy of —0.5 a.u.). Using the LDA for the 
short-range contribution to the exchange energy, rather 
than for the long-range contribution, allows, for a given 
energy, to do a smaller error on this energy, or equiv- 
alently, for a given error, to treat a larger part of the 
exchange energy. Note, in particular, that the LDA er- 
ror on E^'^ becomes vanishingly small toward the left 
end of the plot, corresponding to a very-short-range in- 
teraction. This is in agreement with the fact that the 
LDA becomes exact in the limit of a very short-range 
interaction fll |29| . 

Fig. compares the LDA errors on Ef^, Ef'*, E^ 
and Ept*. It clearly appears that, in the whole energy 
range (except at the end points), the LDA errors on the 
short-range contributions are much smaller (in absolute 
value) than those on the long-range contributions. This 
confirms the appropriateness of the LDA for short-range 
correlations, as often pointed out in the literature (see, 
e.g., Refs. [ll E H3). Note that, contrary to the ex- 
change case, the differences in the relative errors on the 
long-range and short-range contributions are important. 
For example, at an intermediate value of the energy of 
-0.02 a.u., the LDA errors on E 1 *^ and E 1 *^ are -0.0697 
a.u. and —0.0673 a.u., respectively, while the LDA er- 
rors on E s *^ and Ef'^ are as small as —0.0017 a.u. and 
—0.0001 a.u., respectively. One sees in addition that the 
LDA error on E^^ is always significantly smaller (in ab- 
solute value) that the error on Ef'^. We conclude that 
the functional Ef-^ is the best suited for a local density 
approximation. 
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FIG. 2: LDA errors on Ef ' M (thin solid curve), Ef <f (thick 
solid curve), E l £'^ (thin dashed curve) and £^ r ' M (thick dashed 
curve) with respect to -Ef'", -Ef'*, -E]?^ and 
respectively, for the He atom. 

V. CONCLUSIONS 

As regards the long-range/short-range decomposition 
of the exchange functional of Eq. Q), the results of this 
work suggest that the LDA is more accurate for the short- 
range contribution rather than for the long-range contri- 
bution. 



Concerning the long-range/short-range decomposition 
of the correlation functional according to Eq. J2J) or to 
Eq. (JHJ, this work confirms that the LDA is more ac- 
curate for the short-range contributions. The presented 
results suggest in addition that the short-range correla- 
tion functional appearing in the decomposition of Eq. Q 
is better suited for the LDA that the other short-range 
correlation functional appearing in the decomposition of 
Eq. JSJ. 

In the context of the long-range/short-range decom- 
position in DFT, the present paper gives clues about the 
part of the energy, better adapted to (local) density func- 
tional approximations. We hope that it will incite to 
more systematic studies of this topic. 
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